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This paper discusses the progress of work 
which began in mid-2004 sponsored by the 
Office of the Secretary of Defense (OSD) 
Test & Evaluation/Science & Technology 
(T&E/S&T) Program. The purpose of the 
work is to improve the state of the art of 
CFD capabilities for predicting the effects of 
the test media on the flameholding charac- 
teristics in scramjet engines. The program 
has several components including the devel- 
opment of advance algorithms and models for 
simulating engine flowpaths as well as a fun- 
damental experimental and diagnostic devel- 
opment effort to support the formulation and 
validation of the mathematical models. The 
paper will provide details of current work in- 
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volving the development of phenomenological 
models for Reynolds averaged Navier-Stokes 
codes, large-eddy simulation techniques and 
reduced-kinetics models. Experiments that 
will provide data for the modeling efforts will 
also be described, along with with the associ- 
ated nonintrusive diagnostics used to collect 
the data. 

Introduction 

The design and development of a scramjet 
engine is accomplished with several levels of 
analytic tools, ground-based testing and fi- 
nally flight. To achieve the flight conditions 
encountered by a scramjet propelled hyper- 
sonic vehicle in a ground facility, the test gas 
must be heated to high temperatures before 
being introduced into the engine flowpath. 
One method for heating the test gas involves 
mixing fuel, for example hydrogen or butane, 
with the air and allowing it to combust in a 
facility heater prior to use. After the fuel is 
burned, oxygen is added to the test gas to 
make up for the oxygen depleted during com- 
bustion. As a consequence of chemical re- 
action, combustion products, or vitiates, are 
also added to the test gas. These products 
can affect the performance of the engine being 
tested by altering the chemical reactions tak- 
ing place in the engine. Vitiates can change 
the rates of chemical reactions or, in the ex- 
treme, the ability of the reaction to persist 
in the engine combustor. In order to trans- 
late the performance of an engine in a ground 
based facility to the expected performance in 
flight, the effects of facility vitiates must be 
determined. 

This paper will discuss the progress of 
work which began in mid-2004 sponsored 
by the Office of the Secretary of Defense 
(OSD) Test & Evaluation/Science & Tech- 
nology (T&E/S&T) Program. Contributors 
include the NASA Langley Research Cen- 
ter, George Washington University, North 


Carolina State University, the University of 
Pittsburgh, and the University of Virginia. 
The purpose of the work is to improve the 
state of the art of CFD capabilities for pre- 
dicting the effects of the test media on the 
flameholding characteristics in scramjet en- 
gine combustors. The program has several 
components including the development of ad- 
vanced algorithms and models for simulat- 
ing engine flowpaths as well as a fundamen- 
tal experimental and diagnostic development 
effort to support the formulation and valida- 
tion of the mathematical models. The paper 
will provide details of current work involving 
the development of phenomenological models 
for Reynolds averaged Navier-Stokes codes, 
large-eddy simulation techniques, reduced- 
kinetics models, and two experiments with 
the associated nonintrusive diagnostics that 
will provide data for the modeling efforts. 

Axisymmetric Coaxial Free Jet 
Experiment 

Computational fluid dynamics (CFD) 
methods that employ the Reynolds-averaged 
Navier-Stokes (RANS) equations are widely 
used in the design and analysis of hyper- 
sonic airbreathing engine flow paths. These 
methods require models of statistical quan- 
tities of the turbulence in their development 
which have to be empirically calibrated and 
validated. In particular, new models for 
turbulent Schmidt and Prandtl number, as 
well as for turbulence-chemistry interactions, 
are required. 1 While suitable data is avail- 
able for low-speed flows with combustion, 
it is still lacking in supersonic combustion. 
Goyne et al. report measurements using 
particle-imaging velocimetry of mean stream- 
wise velocity in a dual-mode scramjet. 2 At 
the NASA Langley Research Center, several 
data sets have been acquired in a H 2 fu- 
eled supersonic combustor using the coher- 
ent anti-Stokes Raman spectroscopy (CARS) 
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technique 3 * and the dual-Pump CARS tech- 
nique. 4,5 The former technique gave temper- 
ature only while the latter gave both tem- 
perature and composition. Data included 
both mean flow and turbulent statistics, al- 
though the uncertainty in the latter was lim- 
ited by both instrument precision and num- 
ber of measurements from which the statis- 
tics were formed. International work in this 
area includes measurements in scramjet com- 
bustors conducted at ONERA (France) and 
DLR (Germany) using CARS, 6 and other 
non-intrusive techniques. 

Available data are limited to only a subset 
of the important variables (which are temper- 
ature, composition, and velocity) in a limited 
number of geometries, and turbulence data 
are even rarer and of poor precision. To meet 
this need, we are developing an Interferomet- 
ric Rayleigh Scattering technique for measur- 
ing instantaneous velocity to complement our 
Dual-Pump CARS technique. 7,8 This tech- 
nique uses Rayleigh scattering from one of the 
CARS laser beams to measure velocity in the 
same instant that CARS measures tempera- 
ture and composition. Details on the tech- 
niques are given in a later section. Analy- 
sis of streams of such instantaneous measure- 
ments allows us to form the statistical quanti- 
ties (means, variances, covariances) required 
by the modelers. 

Experimental facilities to provide suitable 
flow fields are being developed. An axisymet- 
ric coaxial free jet has been selected since it 
provides the good optical access required for 
the Rayleigh technique, and symmetry can be 
taken advantage of to minimize the number 
of spatial points required. (In order to form 
accurate statistics we need large numbers of 
measurements at each given point - the use 
of symmetry allows the total number of mea- 
surements to be made manageable.) Another 
requirement of the flow field is that it should 
provide data that is relevant to both H 2 - and 


hydrocarbon-fueled hypersonic air-breathing 
engines and the testing of such engines in 
ground facilities that employ vitiated air (the 
products of combustion of either H 2 or a hy- 
drocarbon, enriched with 0 2 to the same con- 
tent as in air). The facility should be capable 
of providing flows of various Mach numbers, 
including supersonic, to establish the effects 
of compressibility. Finally, it is desirable to 
be able to create both supersonic combusting 
flows in which the flame is attached to the 
burner (flame held) and flows in which the 
flame is detached, since both types of flow are 
observed in hypersonic engine combustors. 

Two sets of experiments are being per- 
formed. The second have essentially the same 
flow fields as the first, but are 6.35 times 
larger in scale. The first are conducted in the 
combustion laboratory in which the CARS 
and Rayleigh scattering techniques are being 
developed. The second will be conducted in 
NASA Langley’s Direct Connect Supersonic 
Combustion Test Facility, and will provide 
data in which the turbulent eddy structures 
are better resolved spatially. In this paper we 
will discuss the experiments being conducted 
in the combustion laboratory. 

Figure 1(a) shows the burner and noz- 
zle, sectioned along the axis, without bolts, 
tubes, spark plug and other fittings shown; 
Figure 1(b) is an image of the facility near 
the nozzle exit during operation, showing the 
supersonic jet of vitiated air and the laser 
beams of the CARS-Rayleigh system. The 
facility consists of a water-cooled combustion 
chamber (burner), a silicon carbide nozzle 
(sonic convergent, or supersonic convergent- 
divergent for M=1.6 and M=2), with an exit 
diameter of 10 mm and a coflow nozzle. Var- 
ious combinations of gas flows to the burner 
are possible. In one set of experiments, H 2 or 
CH4 fuel, air and 0 2 are reacted to provide 
vitiated air at various temperatures (depen- 
dent on flow rates). In these experiments, 
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(a) 



Figure 1: The supersonic jet combustor and nozzle, (a) Sectional view (without bolts, tubes, 
etc.), (b) Image during CARS-Rayleigh optical system data acquisition. 


the coflow is of unheated H 2 or CH±, and 
the result (depending on temperature) is a 
mixing and reacting coaxial jet flow. In an- 
other set of experiments, the gas flows to the 
burner are H 2l 0 2 or N 2l and air in such a ra- 
tio that the products contain excess H 2 . The 
coflow is of air and the result again is a mix- 
ing and reacting coaxial jet flow. Depending 
on temperature and on the co-flow rate, the 
flame may be held at the annular base region 
formed between the central and coflow nozzle 
exits where the flow recirculates. 

Computational fluid dynamics (CFD) was 
used to help design both experiments. The 
first step of this effort was for the CFD spe- 
cialist and experimentalist to work together 
to define a set of guidelines to ensure that the 
experiment provided a complete set of data 
that was useful for CFD model validation. 9 
(Often times experimental data sets are in- 


complete, and as a result are of limited value 
for CFD validation.) This effort required the 
experimentalist to identify which properties 
could be measured, where they could be mea- 
sured, and provide estimates of measurement 
accuracy. Similarly the CFD specialist iden- 
tified the properties that were needed, where 
they were needed, and estimates of required 
accuracy. The CFD requirements were based 
on the CFD models to be validated in this ef- 
fort (the turbulent mixing of chemical species 
and thermal energy). The experiment was 
then designed to capture the relevant physics 
of the models to be validated, to limit the 
complexity of the flowfielcl in order to exclude 
physical processes not relevant to the models 
being validated and to limit the sensitivity of 
the flowfielcl to properties which could not be 
readily measured. 

As part of the design effort, a paramet- 
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(a) No Co-flow Injection 
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(b) 15° Co-flow injection 

Figure 2: OH contours showing the sensitiv- 
ity of the flame location to the presence of 
angled co-flow injection. 


ric study of the proposed geometry and run 
conditions was conducted using CFD to de- 
termine the sensitivity of important experi- 
mental quantities to a number of variables in- 
cluding: inner jet temperature, flame holding 
base size, co-axial jet injection angle, inner 
jet Mach number and sensitivity to compu- 
tational values of the turbulent Prandtl and 
Schmidt numbers. An example calculation 
using the Vulcan CFD code 10 is shown in Fig- 
ure 2, which shows the sensitivity of the flame 
location to the presence of a co-flow. This 
flow is with a Mach 2, hot hydrogen-rich in- 
ner jet without and with cold air co-flow. 

Once the experimental geometry was final- 
ized, CFD was used to help define the exper- 
imental test matrix. This matrix, given in 
Table 1, includes variations in the inner jet 
Mach number, the temperature of the inner 
jet, the type of fuel (hydrogen or methane), 
the equivalence ratio of the fuel and air and 
the location of the fuel and air streams (in- 
ner jet or co- flow jet). The CFD calculations 
completed in this step will provide a base- 
line, using the standard constant Prandtl and 
Schmidt number models, to compare with the 
variable Prandtl and Schmidt number models 
described later in this paper. 

Figure 3 contains tables and flow images 
showing the various types of flames that are 
observed in the actual experiment. Two types 
of images are shown: visible light (true color 
images) with a long exposure time, and false 
color infrared (IR) images acquired in the 8 
- 9 micron (long wave) region at an exposure 
time of 10 ms. The tables contain informa- 
tion pertaining to the state of the flame: no 
flame, detached flame, flame held at the base 
or at the external coflow boundary, and if the 
flame holding is marginal, i.e., at the point 
of extinction. Figure 3(a) pertains to cases 
with // 2 -vitiated air in the center jet and H 2 
in the coflow. The center jet sensible enthalpy 
is varied by adjusting the flow rates from that 
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Inner Jet 

Outer Jet 

Mach No. 

Heater Operation 

Mach No. 

Gas 

1,2 

H2+O2+A1T vitiated no 

unreacted H 2 

one T 0 

off 

- 

1,2 

H2+O2+A11 vitiated 

II> rich 

various T a 

off 

- 

1,2 

H2+O2+A1T vitiated 

II) rich 

various T 0 

<1 

Air 

1,2 

H2+O2+A1T vitiated 

O 2 rich 

various T 0 

<1 

ch 4 

1,2 

CH4+O2+A1T vitiated 

O 2 rich 

various T 0 

<1 

ch 4 

1,2 

H2+O2+A1T vitiated 

O 2 rich 

various T 0 

<1 

Hi 

1,2 

CH4+O2+A1T vitiated 

O 2 rich 

various T 0 

<1 

h 2 


Table 1: The experimental test matrix 



Mexlt = 1 

M exit _ 2 

enthalpy M 

$ = 0.5 
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II 

0 = 0.5 

^ = 1.0 

5.5 


EAF 
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EAF 

N 

^ M-DF 

6.5 


EAF 

DP'" 

DF 

7 


EAF 

' DF 

DF 



M exit ~ 1 

M exit 2 

T exit exit -MM exjt -2) 

33% h 2 

50% H 2 

33% h 2 

50% H 2 

1400 K/900 K 


BAF 



1500 K/1000 K 
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BAF 
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M-DF 

1600 K/1100K 
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N = no flame, DF = detached flame, EAF = external attached flame, M = marginal N = no flame, DF = detached flame, BAF = b 

(a) (b) 


attached flame, M = marginal 


Figure 3: Flame state for a matrix of supersonic operating conditions at sonic and Mach 
2 pressure matched exit conditions. Mixed infrared (false color) and visible light images 
(two right-most images), (a) Vitiated air center jet with subsonic hydrogen air coflow at an 
overall equivalence ratio (f). ( b) Hydrogen rich center jet with sonic coflow of air. 
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of Mach 5.5 flight to that of Mach 7, and the 
exit Mach number of the center jet is either 
1.0 or 2.0. The coflow is subsonic at the noz- 
zle exit and the equivalence ratio (the ratio of 
H ‘2 coflow rate to that required to consume all 
the O 2 in the center jet, <p) is either 0.5 or 1. 
(It is not implied that the coflow H 2 reacts 
only with the 0 2 in the center jet.) For M exit 
= 1, the flame is held at the exterior bound- 
ary in all cases, whereas, for M exit = 2, in 
some cases there is no flame and in others the 
flame is detached (stands off from the nozzle 
exit), but there are no cases with flame hold- 
ing. With a detached flame, the trend with 
increasing center jet enthalpy is for the flame 
to move towards the nozzle exit. Figure 3(b) 
pertains to cases with excess H 2 in the center 
jet, either 33% or 50% by volume of the jet 
flow being H 2 . The exit static temperature is 
increased in increments of 100 K, and the exit 
Mach number is either 1.0 or 2.0. The coflow 
of air is sonic and pressure matched at the 
exit. In some cases, no flame is observed; in 
others, detached flames or flames attached at 
the base. The trend with increasing T ex n and 
decreasing M exit is from no flame to detached 
flame to base-attached flame. 

Simultaneous CARS and 
Interferometric Rayleigh 
Scattering 

We have recently reported 11 and summa- 
rize here, the combination of a dual-pump 
coherent anti-Stokes Raman scattering sys- 
tem with an interferometric Rayleigh scat- 
tering system (CARS-1RS) to provide time- 
resolved simultaneous measurement of mul- 
tiple properties in combustion flows. Time- 
resolved simultaneous measurement of tem- 
perature, absolute mole fractions of N 2 , 0 2 , 
and H 2 , and two components of velocity in 
a Hencken burner flame were performed to 
demonstrate the technique. 


The experimental arrangement of the com- 
bined system is shown in Fig. 4. For the 
measurement of temperature and the ab- 
solute mole fractions of N 2 , 0 2 , and H 2 we 
use a dual-pump CARS method. 4 The sys- 
tem uses spectrally-narrow green (injection 
seeded Nd:YAG at 532 nm) and yellow (NB 
Dye laser at 552.9 nm) laser pump beams and 
one spectrally broad red laser (BB Dye laser 
at 607 nm) beam as the Stokes beam. 

The beams are combined at the focusing 
point of a spherical lens Lc (focal length of 
410 mm) in a folded BOXCARS geometry to 
probe Raman transitions of N 2 , 0 2 and H 2 . 
The input beams plus the coherent blue sig- 
nal beam at 491 nm are collected and col- 
limated by another spherical lens with the 
same focal length as Lc. The three input 
beams are captured in a beam dump while 
the blue signal beam is passed to a spectrom- 
eter. The CARS signal, which is a spectrally 
broad blue beam that contains N 2 , 0 2 and H 2 
spectra, is analyzed by a spectrometer and 
recorded by the CCD1 camera. The shape 
of these spectra provides information on the 
temperature while the relative intensities of 
these spectra provide a measure of the rela- 
tive mole fractions. The spectrum is fit with 
a theoretical model to determine the temper- 
ature and mole fractions. 

The velocity measurement is performed 
simultaneously using an interferometric 
Rayleigh scattering measurement system. 7 
The same pulsed, seeded green laser beam 
employed in the CARS technique is used as 
a narrow-band light source for the Rayleigh 
scattering system. The receiving optics 
for Rayleigh scattering are designed to 
capture the Rayleigh scattered light from 
the green pump beam in the measurement 
volume while preserving the scattering angle 
information, and to mix it together with 
the unshifted light of the laser before it is 
passed through the interferometer. In this 


7 

2006 Annual 1TEA Technology Review 



Technology for High Speed and Hypersonic T&E 



Seeded - (tcj:YAG laser @ 20 Hz 
k = 532 nm 


BB Dye laser 
k = 607 nm 


Dove prism 


OF(from reference laser) 


- Perot 


CCD 2 
(Rayleigh) 


Beam stop 

McPherson 

Spectrometer 




NBF 


CARS signal 


CCD 1 (CARS) 


Figure 4: Experimental setup of the combined CARS - IRS system 


way more than one component of velocity 
at multiple points can be measured in one 
interferogram.'’ 12 The velocity components 
measured are those that bisect the angles of 
laser’s incidence and the collection optics. 
Since two directions are collected, two veloc- 
ity components are measured. The Rayleigh 
scattered light from the two measurement 
directions are combined with a polarizing 
beam splitter (PBS) and directed through a 
Fabry-Perot etalon and onto a CCD camera 
for analysis (CCD2). The particular setup 
used here gives a range of measurable veloci- 
ties from 0 to ~3 km/sec up to temperatures 
of about 2500 K. The CARS spectra and the 
Rayleigh spectra are acquired simultaneously 
by synchronizing the cameras with the green 
laser Q-switch pulse at 20 Hz. The spectra 
are subsequently processed, as described in 
reference 11 to determine the temperature, 
composition and velocity. 


To demonstrate the method, experiments 
were carried out in a Hencken burner, which 
provides an adiabatic H 2 -air atmospheric 
pressure flame. The flame was stabilized with 
a co-flow of N 2 . This flame was used because 
it provided a challenging high-temperature 
measurement environment while producing a 
known (near-zero) velocity. 

Simultaneous CARS and Rayleigh scatter- 
ing spectra up to 1610 K are shown in Fig. 5 
(a-c) and in Fig. 5 (d-f), respectively: (a) 
and (d) are in the co-flow of N 2 , (b) and (e) 
are in a high temperature flow containing N 2 
and 0 2l and (c) and (f) are in a high tempera- 
ture flow containing N 2 and a low proportion 
of 0 2 . The CARS spectral plots show the 
experimental data, the theoretical fits, and 
also the residual between them. These spec- 
tra were used to calculate the vibrational gas 
temperature and the mole fraction of H 2 , 0 2 , 
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Figure 5: Simultaneous spectra of CARS (left) and Rayleigh scattering from two viewing 
directions (right). The CARS spectra are both data and fits of theory to the data. The 
Rayleigh spectra are data points connected by lines 


and N 2 . 

For the Fabry-Perot spectra shown in 
Fig. 5, only one fringe order from each spec- 
trum has been analyzed for this paper. The 
narrow shaped peak, more visible in the 
spectra of higher temperature, is the refer- 
ence laser frequency (no Doppler shift). The 
broader spectrum is the Rayleigh scattered 
light. Each figure contains two Rayleigh and 
two reference spectra corresponding to the 
two collecting directions VI (black trace) and 
V2 (red trace of smaller amplitude), differing 
by 34 degrees. Two-component velocity mea- 
surements in the range of 300-1600 K showed 
near-zero velocities (< 30 m/s) and standard 
deviations of 30-40 m/s. These errors are 
about 1% of the dynamic range of this mea- 
surement system (3,000 m/s). Measurement 
of one velocity component was shown to be 
possible at up to ~2,400 K with this system. 
These accuracies and precisions are within 


the desired range required for the planned 
supersonic combustion experiments where ve- 
locities will be up to 1,500 m/s. However, it is 
anticipated that future improvements in both 
hardware and software will allow these errors 
to be reduced by a factor of two or more. 

A Response Surface 
Methodology Approach to 
Refine Computational Models 
Using Experimental Data 

Improved numerical simulation of hyper- 
sonic engine performance relies on the de- 
velopment of enhanced codes with an in- 
creased capability to model turbulence, tur- 
bulent mixing, and kinetics. The validation 
phase of these numerical models requires the 
comparison of simulated results to experi- 
mental data obtained from multiple experi- 
mental cases. One aspect of validation, re- 
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ferred to as model refinement, involves the se- 
lection of model tuning parameters to achieve 
the best agreement with the experimental 
data. 13 

In this section, we propose a response sur- 
face methodology (RSM) approach to select 
the values of the model tuning parameters. 
For over 50 years, 14 RSM has been used to 
successfully characterize and optimize indus- 
trial and scientific products and processes. 
In general, RSM is a strategy to efficiently 
identify important factors in an experimental 
system and to specify the settings of those 
factors that improves a measure of system 
performance. It is applicable to both phys- 
ical experimental systems and computational 
simulation. 

The model refinement phase is essentially 
an experiment in which we set specific values 
of the model tuning parameters, execute the 
simulation, measure the correlation to exper- 
imental results, and seek to find the values of 
the parameters that improve the correlation, 
or agreement, with actual experimental mea- 
surements. In the RSM context, the model 
tuning parameters are the factors with as- 
sociated levels, which represent their numer- 
ical values. The simulation-to-experimental 
correlation forms the response that we seek 
to improve. Applying RSM offers a gen- 
eral, systematic approach to perform model 
refinement that emphasizes the use of mini- 
mal computational resources and features an 
analysis approach to gain deeper insights into 
the underlying physics. For more information 
on statistical experimental design and RSM, 
see Refs. 15,16. 

The first step in the RSM approach is to 
choose the model tuning parameters (factors) 
to be varied and their high and low levels, or 
range. Then, an experimental design is con- 
structed that specifies the levels and combi- 
nations of factors to be run through the sim- 


Model Tuning Parameters 


Run 

Cy 

Cy-2 

7 

CO 

C'Y4 1 

Cyj 

1 

-1 

-1 

-1 

-1 

1 

2 

1 

-1 

-1 

-1 

-1 

3 

-1 

1 

-1 

-1 

-1 

4 

1 

1 

-1 

-1 

1 

5 

-1 

-1 

1 

-1 

-1 

6 

1 

-1 

1 

-1 

1 

7 

-1 

1 

1 

-1 

1 

8 

1 

1 

1 

-1 

-1 

9 

-1 

-1 

-1 

1 

-1 

10 

1 

-1 

-1 

1 

1 

11 

-1 

1 

-1 

1 

1 

12 

1 

1 

-1 

1 

-1 

13 

-1 

-1 

1 

1 

1 

14 

1 

-1 

1 

1 

-1 

15 

-1 

1 

1 

1 

-1 

16 

1 

1 

1 

1 

1 

17 

0 

0 

0 

0 

0 


Table 2: Initial Experimental Design (+/-1 
represent high and low levels of each factor 
in coded units, 0 represents the mid-point) 

illation model. For each combination of fac- 
tors, known as an experimental run , measures 
of simulation-to-experimental agreement (re- 
sponses) are obtained. From these data, 
parametric validation models of the relation- 
ship between the factors and the responses 
are estimated. These validation models de- 
scribe a multidimensional validation surface 
of the difference between the simulated and 
experimental results as a function of the 
model tuning parameters. Our goal is to 
identify the regions of the validation surface 
(combination of factors) that represent agree- 
ment with the experimental results. Since the 
functional form of the validation surface is 
unknown, an iterative process is performed 
to assess the adequacy of an estimated val- 
idation model and augment the experimen- 
tal design as required, thereby enabling the 
estimation of a higher fidelity model. Once 
adequate models are found, they are com- 
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bined to perform multiple response optimiza- 
tion, thereby enabling the determination of 
the values of the model tuning parameters to 
achieve the best correlation to the experimen- 
tal data for all of the response quantities of 
interest. Note that we use the term best to 
describe a trade-off among multiple compet- 
ing criteria in correlating different aspects of 
the simulation model. 

In the final step, the values of the model 
tuning parameters determined by the opti- 
mization are run through the simulation code 
to confirm that the predicted quality of cor- 
relation with the experimental results is ob- 
tained. This confirmation phase provides 
confidence in the estimated validation models 
and the optimization results. 

As an example of the initial setup of this 
proposed approach, consider the development 
of a variable turbulent Schmidt number for- 
mulation, 1 ' described in the next section, to 
model a Helium co-flow mixing experiment. 18 
Five model tuning parameters ( Cy , CQ 2 , CVs, 
Cy 4 i , Cyt) that proportionally scale specific 
terms within the turbulence equation, are 
varied at three levels; a high, middle, and 
low setting. The initial experimental design 
requires 2 5_1 = 16 factorial combinations of 
the five factors at their high and low levels, 
representing a one-half fraction of all possi- 
ble combinations. In addition, the design in- 
cludes one combination with all 5 factors si- 
multaneously set at their middle level. The 
design is provided in Table 2 in coded units. 
Geometrically, the design is a subset of ver- 
tices of a five-dimensional hypercube and a 
design point in the center of the hypercube. 

This design supports the efficient and 
unique estimation of a first-order model given 
by 

k k — 1 k 

y = Po + 52 P iXi + 52 52 Pij x i x j + e, 

i = 1 i = 1 j=i -\- 1 

where y is the response (validation metric), x 


is the level of the factors (model tuning pa- 
rameters), /3’s are the validation model coef- 
ficients, e is the residual error in fitting the 
model, k is the number of parameters (fac- 
tors) in the experiment. Typically in these 
types of simulation codes, £ is a systematic 
error rather than a random error clue to the 
deterministic nature of the simulation. 

An important feature of this model is that 
it allows for the modeling of interactions be- 
tween the factors, represented by the 
terms. If the effect on the response by chang- 
ing the value of one of the model tuning pa- 
rameters depends on the setting of another 
parameter value, then we say that the factors 
interact. Synergistic and antagonistic inter- 
action effects are common in many applica- 
tions and cannot be captured if the factors 
are varied one at a time rather than simul- 
taneously in a factorial arrangement. Identi- 
fying interactions often provides new insights 
into the underlying physics of the experimen- 
tal system. 

In this example, the four system response 
quantities 19 (SRQ) considered are species 
mass fraction, velocity, Pitot pressure, and 
total temperature measured throughout the 
flow field. Typically, the correlation to exper- 
imental results involves a subjective, qualita- 
tive review of plots. However, a quantitative 
measure of correlation quality is required that 
has the form of VM = / (y c — y e ), where 
VM is the validation metric vector, y c is the 
SRQ from the computational model, and y e 
is the SRQ from the experimental data. Note 
that y is a vector containing the values of 
the responses at the location in the flow field 
where experimental measurements were ob- 
tained. At this stage in the development of 
the approach, a simple L 2 norm of VM is 
proposed to provide a numerical measure of 
the agreement with the experimental results. 
A full development of this case study will be 
presented in a subsequent paper. 
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In summary, the proposed RSM approach 
to model refinement is expected to offer a gen- 
eral, structured approach to obtain values for 
the model tuning parameters. In addition, 
due to the structured nature of the approach, 
the selection of the tuning parameters will be 
reproducible by other researchers. A particu- 
lar strength of the approach is its straightfor- 
ward extension to higher- dimensional factor 
spaces and its ability to incorporate multiple 
experimental cases, thereby providing a set of 
parameters that are adequate over a range of 
experimental conditions. 

Improved Turbulence Models 
for Reacting Flows 

The research to improve the modeling of 
turbulent transport of heat and mass has re- 
sulted in a tentative turbulence model where 
the turbulent Prandtl (Prt) and Schmidt 
(Set) numbers are calculated as part of the 
solution and where reactions involving chem- 
ical source terms are modeled. The turbu- 
lent Prandtl and Schmidt numbers are ob- 
tained by solving additional equations for the 
variance of enthalpy and its dissipation rate; 
and for the variance of concentrations and 
their dissipation rate. The resulting equa- 
tions, together with modeling of terms involv- 
ing chemical source terms, are given in the 
references . 17 It is to be noted that the model 
avoids the use of assumed or evolution PDF’s 
to account for chemistry/turbulence interac- 
tions. 

Two sets of H 2 /Air chemical kinetic 
mechanisms are considered: the seven 

species/seven reaction model of Jachimowski 
, 21 where the reaction rates are functions of 
temperature; and the nine species/nineteen 
reaction model of Connaire, et al. 22 where 
the reaction rates are functions of both pres- 
sure and temperature. 

The predictions from the model are com- 
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H2 Injection 


Figure 6: Schematic of Experiment Setup 


pared with two sets of experiments involv- 
ing supersonic mixing 23,24 and combustion. 24 
Only results involving combustion will be pre- 
sented here. Results involving mixing are 
given in Ref. 17. A schematic of the experi- 
ment of Burrows and Kurkov 24 is shown in 
Fig. 6. Hydrogen is injected into the test sec- 
tion at Mach M=1 through a nickel injector 
parallel to the vitiated main flow. At the 
entrance of the test section, M = 2.44, the 
state pressure is one atmosphere, and the sta- 
tic temperature is in the range 1250-1270K. 
Hydrogen is injected at a total temperature 
slightly above the ambient temperature. All 
comparisons are made at the test section exit 
plane, which is 35.6 cm from the hydrogen in- 
jection step. Three sets of figures are pre- 
sented. For the first two sets, Jachimowski’s 
model is employed. In the first set, terms 
involving averages of chemical source terms, 
which represent the turbulence-chemistry in- 
teractions, are ignored, while in the second, 
the contributions of these terms are included. 
As is seen from Fig. 7, poor agreement with 
experiment is indicated. Figures 8 and 9 show 
contours of Sc t and Pr t . It appears that the 
main source of the discrepancy is a result of 
reduced Pr t near the mixing region. This 
has the tendency of promoting heat trans- 
fer and early combustion. Figure 10 shows 
that when turbulence-chemistry interactions 
are included, much better agreement with ex- 
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Pr . 



Figure 9: Prandtl Number Contours, with- 
out Chemical Source Term, Jachimowski’s 7 
Species/7 Reactions Model 


Figure 7: Comparison of Computed and Mea- 
sured Volume Fraction, with Chemical Source 
Term, Jachimowski’s 7 Species/7 Reactions 
Model 
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Figure 8: Schmidt Number Contours, with- 
out Chemical Source Term, Jachimowski’s 7 
Species/7 Reactions Model 



Figure 10: Comparison of Computed and 

Measured Volume Fraction, with Chemical 
Source Term, Jachimowski’s 7 Species/7 Re- 
actions Model 
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Figure 11: Schmidt Number Contours, with 
Chemical Source Term, Jachimowski’s 7 
Species/7 Reactions Model 


Figure 13: Comparison of Computed and 

Measured Volume Fraction, with Chemical 
Source Term, Connaire et al 9 Species/19 Re- 
actions Model 
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Figure 12: Prandtl Number Contours, with 
Chemical Source Term, Jachimowski’s 7 
Species/7 Reactions Model 


periment is indicated. Figures 11 and 12 show 
contours of Sc t and Pr t . In the third set, 
Figures 13-15 show the prediction of Con- 
naire 22 et al. model. As is seen from Fig. 13, 
good agreement with experiment is indicated. 
This suggest that modeling approaches can 
serve as a substitute for approaches employ- 
ing evolution PDF’s. The differences in the 
Sc t contours of Figs. 11 and 14 illustrate the 
fact that, in regions where mixing does not 
play a role, the value of diffusion coefficient 
has little influence on the final results. 

Based on the above results, two relevant 
observations can be made. The first is 
that modeling of averages of terms involv- 
ing chemical source terms is a viable op- 
tion. When this approach is compared with 
approaches requiring assumed or evolution 
PDF’s, a great deal of computational effi- 
ciency is achieved. Second, a relatively in- 
expensive calculation of variable Pr t and Set 
can be obtained by assuming a value for 
the Lewis number, and eliminating either the 
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Figure 14: Schmidt Number Contours, with 
Chemical Source Term, Connaire et al 9 
Species/19 Reactions Model 
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Figure 15: Prandtl Number Contours, with 
Chemical Source Term, Connaire et al 9 
Species/19 Reactions Model 


equation for the enthalpy variance and its 
dissipation rate, or those for the variance of 
concentrations and its dissipation rate. This, 
however, will result in ignoring one of the av- 
erages involving chemical source terms. Be- 
cause inclusion of such terms is important, 
assuming a constant Lewis number is not rec- 
ommended. 

Large Eddy Simulation of 
Reacting Flows 

It is now widely recognized that one of the 
most convenient means of predicting the un- 
steady physics of turbulent reacting flows is 
via large eddy simulation (LES). 34 A serious 
issue in such simulations is accurate mod- 
eling of the subgrid scale (SGS) quantities. 
While this issue is important in any LES, it 
is particularly difficult when dealing with re- 
acting flows . 28 The filtered density function 
(FDF) has proven particularly effective for 
accurate SGS modeling in reacting flows. 32 
This method is the counterpart of the prob- 
ability density function (PDF) methodology, 
which has been long-established in RANS cal- 
culations of turbulent combustion , 38 A re- 
liable means of determining the FDF is via 
solution of its transport equation. A system- 
atic study in this endeavor was initiated by 
Colucci et al. 29 who developed and solved a 
modelled transport equation for the marginal 
scalar FDF (SFDF) in constant density re- 
acting flows. The extension of this work was 
conducted by Jaberi et al. 33 who developed 
the marginal scalar filtered mass density func- 
tion (SFMDF). This is essentially the mass- 
weighted form of the SFDF. Due to the mar- 
ginal nature of the FDFs in these studies, all 
of the hydrodynamic effects including all of 
the velocity-scalar correlations are modelled 
via conventional (non-FDF) methods. The 
extension for modeling of the velocity field 
was done by Gicquel et al. 31 who developed 
and solved a transport equation for the mar- 
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ginal velocity FDF (VFDF) in constant den- 
sity flows. This work demonstrated some of 
the advantages of the FDF in comparison to 
conventional methods in accounting for the 
effects of SGS velocity correlations. The most 
general formulation, is the joint FDF of the 
velocity and the scalar fields (VSFDF) in con- 
stant density flows. 41 

The objective of this part of the work is to 
extend the capabilities of the FDF method- 
ology for LES of the flows considered in the 
other parts of this investigation. To achieve 
this objective, our work has been concen- 
trated on two main tasks: (1) Development of 
the joint velocity-scalar filtered mass density 
function (VSFMDF), and (2) Application of 
the SFMDF for prediction of complex turbu- 
lent combustion systems. The work in (1) 
was motivated simply because of the physics 
of high speed turbulent combustion which re- 
quire variable density formulation. The work 
in (2) was motivated in our continuing at- 
tempt to extend the boundaries of applica- 
bility of the FDF method. 

In the work pertaining to (1), an ex- 
act transport equation is derived for the 
VSFMDF. It is shown that the effects of SGS 
convection and chemical reaction appear in 
closed form in this equation. The unclosed 
terms in this equation are modelled in a fash- 
ion similar to that typically used in RANS. A 
set of stochastic differential equations (SDEs) 
are considered which yield statistically equiv- 
alent results to the modelled VSFMDF trans- 
port equation. The SDEs are solved numeri- 
cally by a Lagrangian Monte Carlo procedure 
in which the Ito-Gikhman character of the 
SDEs is preserved. The consistency of the 
proposed SDEs and the convergence of the 
Monte Carlo solution are assessed. In non- 
reacting flows, it is shown that the VSFMDF 
results agree well with those obtained by 
a “conventional” finite-difference LES proce- 
dure in which the transport equations corre- 


sponding to the filtered quantities are solved 
directly. The VSFMDF results are also com- 
pared with those obtained by the Smagorin- 
sky closure, and all the results are assessed 
via comparison with data obtained by direct 
numerical simulation (DNS) of a temporally 
developing mixing layer involving transport 
of a passive scalar. It is shown that all the 
first two moments including the scalar fluxes 
are predicted well by the VSFMDF. In these 
simulations, the VSFMDF methodology is 
shown to be able to represent the variable 
density effects very well. The VSFMDF pre- 
dictions of shear layer growth rate for various 
density ratios compare well with DNS data. 
The predictive capabilities of the VSFMDF 
in reacting flows are further demonstrated by 
LES of reacting shear flows. The predictions 
show favorable agreements with laboratory 
data, 36 and demonstrate several of the fea- 
tures as observed experimentally. For details 
of this work, we refer to a recent Ph.D. Dis- 
sertation by Sheikhi 39 

In the work pertaining to (2) the SFMDF 
is applied for LES of a turbulent bluff-body 
stabilized hydrogen- methane jet flame. This 
flame is studied by the Combustion Research 
Facility at the Sandia National Laboratories 
and the Thermal Research Group at the Uni- 
versity of Sydney. 25, 26,35 This work was pro- 
moted in the wake of our recent success 40 
in SFMDF predictions of the simpler piloted 
jet flame experiments of Sandia . 25,27,3 ‘ For 
this flame (which exhibits little local extinc- 
tion), a simple flamelet model is used to relate 
the instantaneous composition to the mix- 
ture fraction. The modelled SFMDF trans- 
port equation is solved by a hybrid finite- 
difference/ Monte Carlo scheme. The results 
via this method capture important features 
of the flame as observed experimentally. In 
fact, all of the results including the means, 
RMS values and PDFs of the scalar field as 
predicted by SFMDF compare very well with 
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experimental data. For details of this work, 
we refer to a recent Ph.D. Dissertation by 
Drozda 30 

Reduced Chemical Kinetics 
Models 

While turbulent transport and turbulent 
reaction models are system dependent (eg. 
inflow turbulent intensity), the fundamental 
laminar reaction models can be investigated 
independent of the actual reaction flow field 
by analyzing the models under relevant phys- 
ical time scales. In particular, the present ef- 
fort is focused on the development of reduced 
reaction models based on laminar kinetic ef- 
fects on ignition/propagation/extinction of 
hydrocarbon-air mixtures and subsequent ap- 
plication to turbulence-chemistry models. 

Because of the importance of ethylene in 
many hydrocarbon oxidation pathways, both 
elementary 42 as well as ad hoc 4 3 models have 
been developed to describe the ignition and 
oxidation processes. Based on the elementary 
models, systematically derived reduced reac- 
tion models have also been developed for en- 
gineering applications; however, these models 
have focused on describing either the propa- 
gation of premixed flames and extinction of 
nonpremixed flames 44 or ignition. 45 In gen- 
eral, reduced reaction models developed for 
flame propagation and extinction cannot be 
applied to the prediction of ignition phenom- 
ena because some of the key initiation ele- 
mentary steps necessary for ignition phenom- 
ena may be missing or transport effects are 
not included in ignition analysis. As part of 
the present work, a general methodology has 
been developed to address above shortcom- 
ings, as described below. 

In early investigations using reduced reac- 
tion models based on the quasi-steady-state 
(QSS) approximation, the computational ef- 
ficiency was not a primary factor in devel- 


oping these reduced models. For example, 
highly-stiff coupled algebraic relations arising 
from the QSS approximation were solved it- 
eratively, i.e. with “inner” iterations, or the 
coupled QSS algebraic relations were trun- 
cated to yield “explicit” expressions for the 
species in QSS, albeit within a narrow range 
of validity or with diminished accuracy. 46 
The formal computer based reduction meth- 
ods, e.g. computational singular perturba- 
tion (CSP) method of Lam, 47 computer-aided 
reduction method (CARM) of Chen, 48 and 
intrinsic low-dimensional manifold (ILDM) 
method of Pope, 49 facilitated the develop- 
ment of reduced reaction models for large 
detailed models. However, these mathemat- 
ical/computational approaches still require 
the solution of highly-stiff algebraic relations 
using computationally expensive “inner” it- 
erations. The automatic mechanism reduc- 
tion procedure developed here (based on the 
MATLAB programming language), can re- 
lax certain QSS relationships to obtain “ex- 
plicit” expressions for the QSS species, i.e. 
avoid costly “inner” iterations. The level 
of relaxation can be based on a predeter- 
mined threshold satisfying the QSS approx- 
imation. LInlike other automated reduction 
procedures, the present explicit reduced reac- 
tion modeling (ERRM) approach can be ap- 
plied to ignition phenomenon as well as pre- 
mixed and non-premixed flames . 50 In this re- 
gard, the present MATLAB-based ERRM ap- 
proach is more general compared to the CSP 
approach. 47 

Figures 16 and 17 show a typical compar- 
ison of the ignition delay for C 2 H 4 - air and 
H 2 - air mixtures, while Fig. 18 shows a com- 
parison of premixed flame propagation veloc- 
ity predictions for CH 4 - air mixtures based 
on two reduced reaction models developed 
for C 2 H±. Also shown are the results us- 
ing Wang’s 42 detailed kinetic model with 73 
species and 469 reactions and a skeletal model 
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Figure 16: Comparison of ignition delay pre- 
dictions of stoichiometric C2H4-air mixture 
as a function of temperature 

developed consisting of 31 species and 128 re- 
actions. The 15-step reduced reaction model 
is based on flame propagation of C 2 H A - air 
flames, while the 18-step reduced reaction 
model is based on C 2 H 4 - air ignition delay , 50 
To demonstrate the robustness of re- 
duced reaction models, the current 15-step 
and 18-step models are used in predicting 
the flame extinction condition of C 2 H^ and 
air non-premixed flame, with the air-stream 
heated to 900 K using a typical vitiation fa- 
cility (eg. iso-butane preheater with vitiated 
air composition of X 02 = 0.21, X N2 = 0.727, 
Xco 2 = 0.028, and X H2 o = 0.035). As shown 
in Fig. 19, the extinction strain rate is re- 
duced because of the vitiation species C0 2 
and H 2 0 included in the model and the 18- 
step reduced reaction model can capture the 
same trend as the skeletal model. In the fu- 
ture, effects of other intermediate species in 
vitiation facilities on ignition and flame prop- 
agation and extinction will be analysed. 

Concluding Remarks 

This paper has described work involving 
the development of phenomenological models 



Figure 17: Comparison of ignition delay of 
predictions of stoichiometric H2-air mixture 
as a function of temperature 



Figure 18: Comparison of laminar burning 
velocity of premixed CH4-air mixture as a 
function of equivalence ratio, at T 0 =300K 
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Figure 19: Comparison of peak flame tem- 
perature prediction of C2H4 and air non- 
premixed flame with flow strain rate, with air 
preheated to 900 K using a vitiated facility 


for Reynolds averaged Navier-Stokes codes, 
subgrid scale models used in large-eddy sim- 
ulation and reduced-kinetics models to study 
the effects of vitiation on engine testing in 
combustion heated facilities. Two fundamen- 
tal experiments are being performed to pro- 
vide data that will be used in the develop- 
ment and refinement of these models. Ex- 
perimental data is extracted from the experi- 
ments using nonintrusive diagnostics that al- 
low accurate simultaneous measurement of 
temperature, species, and up to three com- 
ponents of velocity in supersonic flow with- 
out changing the character of the flowfield. 
Once the databases are in hand, the data is 
analyzed using a response surface methodol- 
ogy that provides an efficient means of deter- 
mining critical parameters in a chosen model. 
The models will then be incorporated into 
combustion codes used in engine flowpath 
analysis and design. 
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